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Abstract 

Consider a free electron gas in a confining potential and a magnetic field 
in arbitrary dimensions. If this gas is in thermal equilibrium with a reservoir 
at temperature T > 0, one can study its orbital magnetic response (omitting 
the spin). One defines a conveniently "smeared out" magnetization M, and 
the corresponding magnetic susceptibility Xi which will be analyzed from a 
semiclassical point of view, namely when Ti (the Planck constant) is small 
compared to classical actions characterizing the system. Then various regimes 
of temperature T are studied where M and x can be obtained in the form 
of suitable asymptotic ^-expansions. In particular when T is of the order of 
h, oscillations "a la de Haas-van Alphen" appear, that can be linked to the 
classical periodic orbits of the electronic motion. 
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1 Introduction 

The magnetic response theory for a free electron gas is an old problem conside- 
red by Landau [ 15| , Fock and Peierls [20|. The revival of interest in physics arose 



from the advances of recent experiments that made possible measurements of the 
magnetic response on small 2-dimensional electronic devices. These devices are so 
"pure" that the classical as well as quantum motion inside them can be considered 
as "ballistic", i.e. is uniquely determined by the confining potential. (Taking into 
account impurities would consist in adding the random potential created by random 
point scatterers inside the material). The bi-dimensional structure of such electron- 
gas is realized through semi-conductor heterostructures whose size, and shape can 
be controlled experimentally, together with the number N of confined electrons. The 
system being in contact with a reservoir at temperature T, and submitted to a ma- 
gnetic field B perpendicular to the surface, the magnetic response can be measured : 
say, the magnetization M or magnetic susceptibility x as a function of the thermody- 



namic parameters T, N, B. ||22j| . These experiments manifest the sensitivity of the 
magnetic response to the integrable versus non-integrable character of the classical 
dynamics of one electron in the system. Numerical experiments on two-dimensional 
magnetic billards have confirmed this observation, and suggested that the quantum 
magnetic response is an experimentally accessible criterion for distinguishing classi- 
cally integrable versus chaotic dynamics WJ\. A number of theoretical studies have 
analyzed the magnetic response from a "semi-classical" point of view, namely as 
properties manifesting themselves in the limit when h (the Planck constant) is small 
compared to classical actions characterizing the system (say h <C a 2 eB/c where a is 
a typical size of the system, e the charge of the electron, c the velocity of light and 
B the magnetic field size) |£§. In these studies, the Coulomb interactions between 
the electrons in the system are neglected, so that the system is a "free electron gas" 
to which the usual thermodynamic formalism is applied. 

The thermodynamic functions in the grand-canonical ensemble can be expressed 
through the density of states of the quantum Hamiltonian for one electron in the 
system. This quantum density of states, in the semi-classical limit, splits into a 
mean part and a strongly oscillating one, according to the well known semi-classical 
trace formula. This formula is known in mathematics as Poisson formula (Colin 
de Verdiere H , Duistermaat-Guillemin M) and in physics as the Gutzwiller trace 



formula in the chaotic case [12], or the Berry- Tabor trace formula in the integrable 
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case |y. This splitting provides a similar splitting in the magnetic response, which 
allows to understand the oscillations "a la de Haas-van Alphen" of the magnetic 
susceptibility and their link with the classical periodic orbits of the electronic motion. 

The aim of the present paper is to reconsider these questions from a mathematical 
point of view, in the following two directions (for non-interacting electron gases in 
arbitrary dimension, and not necessarily homogeneous magnetic fields) 

- examine the regimes of temperature in which the magnetic response can be 
obtained semi-classically in the form of asymptotic h expansions 

- investigate a "mesoscopic" regime of low temperatures where the periodic orbits 

of the classical one-electron dynamics manifest themselves as highly oscillating 
contributions, to the magnetic response. 



In a recent work, Fournais |TT| studies the semi-classics of the quantum current 
for a non-interacting gas of electrons in dimension n and temperature T, confined 
in a potential V and subject to a suitable magnetic field B. For fixed non-zero 
temperature T, he obtains a complete asymptotic expansion of the quantum current 
in small h, and for zero temperature, he obtains the dominant contribution plus an 
error term under suitable assumptions. J. Butler || has recently reexamined this 
last case using a "semi-classical trace formula" by Petkov and Popov . 



We recall that in all these studies, the spin of the electron is omitted so that 
only the orbital magnetic response is considered. 
The content of our paper is the following : 

• In section 2 we consider the case when the temperature T is large compared 
to the Planck constant h. We prove asymptotic expansion in h for the ther- 
modynamical potential and we recover the Landau diamagnetic formula for 
2-dimensional free electron gas. 



• In section 3 we consider the case where the temperature T is of the same 
order as h. Then we prove that the magnetization splits into two terms : 
an average part with a regular asymptotic expansion in h plus an oscillating 
part in h which is the contribution of the periodic orbits of the classical motion. 
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• In section 4 we come back to the regime T larger than h and prove that the 
contribution of non zero periods of the classical motion is exponentially small 
in h. 

2 The Landau magnetism 

We shall first give the notations and assumptions that will hold all along this 
paper. 

Given (5 > 0, we set : 

F p (x) = ~Log(l + e-^) (2.1) 

fp(x) = F>(x) = (l + eP*)- 1 (2.2) 

fp is related with the Fermi-Dirac distribution. 

These functions are meromorphic, with poles (or cuts for Fp) at : 

2A: + 1 . , _ 
x = — - — nx k G L 

P 

Let k £ IR be a real parameter (coupling constant with a magnetic field). We 
consider a family of Hamiltonians with magnetic fields given by : 

H K (q,p) = ^(p-Ka(q)) 2 + V(q) (2.3) 

where V : IR n i— > IR and a : IR n i— > IR n are C°° functions satisfying the following 
properties. 

(H.l) Vg G IT, \/(g) > 1 , |^\/(g)| < C a y(g) 

(H.2) Vg G R™, |^a(g)| < C a V(g) 1 / 2 

(H.3) Vg G R n , V(q) > c (l + |g| 2 ) S/2 some s, c > 

(confinement assumption). 
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Let now H K be the Weyl quantization of H K . The previous assumptions ensure 
that H K is self-adjoint and its spectrum cr(H K ) C [s, <x>) is pure point for \k\ < n 
where e > (fl23f). Let us call (Ej)j £]N and (</?j)j'eiv the set of corresponding 
eigenvalues and eigenstates. 

In the grand-canonical ensemble, the thermodynamic potential Q is given by : 

n(/3, //, «) = £ F p( E i ~ = 1* (#« " } (2-4) 

where \x > is the chemical potential, /3 = being the Boltzmann 

constant, and T > the temperature, (re will be the size of the magnetic field). 
Furthermore, the mean-number of particles in the grand-canonical ensemble is given 
by : 



N(j3,fjL,K)=Tr[f /3 [H K - f i)) (2.5) 

Using the functional calculus p3| , it is not difficult to see that Fp(H — //) and 
fp(H — /i) are trace-class and that the function : k h- > //, k) is C°° for |/c| < Kq. 
We shall denote d K = -j^. 

Proposition 2.1 T/ie function : k \— > Q(j3,[i, k) is C°° on R. In particular we have 



d K Q = Tr 



(# K - //) (2.6) 



Now we have the following definitions of magnetization M and magnetic suscepti- 
bility x ■ 



M = d K Q 

X 



Tr 

d K M 



fp \H K - fij d K H K 



(2.7) 
(2.8) 



Proof of propositionfO] Since cr(H K ) C [eq, +oo), we can draw a suitable curve A 
in the complex energy plane, around a(H K ), with all branching points of Fp(z — /i) 
left outside. So using Cauchy formula we have 



[z - (jl) [H k - z 



-i 



Using Lebesgue convergence theorem and cyclicity of the trace, we get 



1 

2in 



Tr / dz 



F p {z - n) (H K - z) d K H K 



(2.9) 



(2.10) 
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Integration by parts give 

d K Q = ^Tr jT dzf p (z - (H K - z)' 1 d K H K (2.11) 

This procedure can be easily iterated to prove that Q is C°°-smooth in k. Moreover 
in the semiclassical regime we can prove that the asymptotics for derivatives in k of 
Q can be computed using the following commutators formulas for derivatives of the 
resolvent. Starting from the well known identity : 

[A, (H - zy 1 ] = (H- zy'iH, A](H - z)~\ (2.12) 

we get 

d K (H - z)' 1 = d K H(H - zy 2 - [d K H, H]{H - z)' 3 
-[H, [H, d K H]]{H - z)' A + (H- z)- l [H, [H, [H, d K H]]](H - z)~\ (2.13) 

Each commutator gives one ft and we can compute in the same way higher derivatives 
in k. So using Cauchy formula we can compute asymptotics in ft of derivatives in k 

of a ■ 



Theorem 2.2 For any e > and k,q > 0, f2 admits an asymptotic expansion in ft, 
uniform in k for \k\ < k , and for (3 < % £ ' 2 ^ . More explicitly, for any N E IN we 
have : 

n = h~ n E T , w n jk + o (n N+1 " n (2.i4) 



i=o fc <M 



with 



Q jk = dq dp d jk (q,p) F^\H k - //) 

J M 



djk being a suitable linear combination of derivatives of H K with respect to q, p and 
k(n) a constant depending only on the dimension n (7c (n) <2n+l). In particular : 



^oo= I dq dp Fp (H K (q,p) - fi) 

JR 2n 



(2.15) 



-f2 22 — g^23 



P 



48tt 2 



/ dq dp 

Jm 2 " 



II B(q) \\ 2 -Ed 2 k V 



jk 



cosh 



l(H K (q,p)-t,) 



(2.16) 



where Bjk is the magnetic field 
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B 



j<k 



B 



da* da i 



dlv 



d 2 V 



dq k dqj' 3k 



(2.17) 



i0 dqjdq k 

and we have chosen the gauge so that da/dq is symmetric. Moreover, the asymptotic 
expansion can be derived term by term with respect to k and yields an asymptotic 
expansion of the magnetization and the magnetic susceptibility. 



Proof . We start with the following Cauchy formula as in the proof of proposition 
Ffs -//) = _ jf dz F p {z - y) (H K - z)' 1 



2.18 



Proceeding as in [23|], good enough semi-classical approximations of (H K — z) 1 (for 
z G A) are obtained for any integer iV of the following form : 



-l N -l 
(H K -z) =J2Wb j (z)-h N+1 (H K -z) R N (z) 



(2.19) 



3=0 



where A denotes Weyl quantization, and bj(z) for j > 2 are obtained, from bo(z) 
(H K (q,p) — z)^ 1 by the formula 



(2.20) 



2<£<[f] 



dj£ being a symbol constructed through partial derivatives of H K (q,p) and that can 
be computed explicitly. Furthermore, due to the particular form ( |2.3| ) of H K , bj = 
for odd j's, and Rn obeys : 





z 


\Rn(z)\<C n 


Im z 



(2.21) 



k(n) depending only on the dimension n. (for details see || 
Now inserting ( |2.19 ) into the Cauchy formula and using 



/ (m) (A) 
together with (h = 27th) : 



■l) m m\ 
~2in 



dz f(z) (z - A) 



-m— 1 



Tr A = h' n / dp dq A(q,p) 

J]R 2n 



(2.22) 



(2.23) 



we get the result. Let us make explicit the calculus. We have to find 62(2) such that 
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(H K - z) (b (z) + h 2 h(z)) = 1 + Ofi') (2.24) 
which, according to the rule for the symbol of the product of two operators, yields 



b 2 {z) = d 22 b 3 (z) + d 23 b*(z) 
with b (z) = (H K — z) -1 , and 

d 22 = ~ E (^'H K )(d2'd«H K 

4 \a\+\a'\=2 



(2.25) 



ala ! 



.1)KI 



(2.26) 



|a|+|a'|=2 

a being a multi-index a = (an, ■ ■ ■ a n ) G -CV™, we denote as usually by 9? the multiple 
derivative • • • by |a| the sum a j, an d by a\ the product J] 



J'=l 



Now the calculi proceed as in [13 



i „ / n 

-2 («9 Pfc tf K ) (d 2 p . qk H K ) + (d Pj H K ) (d Pk H K ) (d 2 g . qk H K 

and integrating by parts (over qj, or p^) we get 



Q 23 = 2 / dg dp if } (# K - V)d 22 (q,p) = 2n 22 (2.27) 



and thus : 



O O 



We now make ^22 explicit 



(2.28) 



= J |« 2 II 5 II 2 + E (« (ft - ™M) ■ djM - d k V{qj) I (2.29) 
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Clearly the term (pj — Kdj(q))dj k a(q) does not contribute to ^22 using the change of 
variable p — > p — na(q), and the oddness of the integrand with respect to p variable. 
Thus we are left with 



n 



22 



~ J^Jq dp if> (flo -/*) ^ 2 II f -E^*V(g)j 



Now the uniformity of the asymptotic expansion in h with respect to (3 < ft 
(for any e > 0) comes from the fact that F^\x) = (3 k ~ l F[ k \(3x) , (A; > 1) so that 
fr J '/3 fc < W £ for < Furthermore, the error term in (|2.14j) follows from (2.21). 

For the magnetization M, we start with formula ( j2.7|) and we use the semi- 
classical expansion for the resolvant to get 

M = ^-Tr J dzf (z - y) ( E Wb^dJLA (2.30) 



yO<j<N 



N+l 



Tr / dzf p (z — fi) (H K — z)- 1 R N (z)d K H K (2.31) 
Am J A 

Then using integration by parts we can prove that the first term in ( |2.30| ) is 



h ~ n E E Mr-^A*' ( 2 - 32 ) 

0<j<Nk<3j/2 

and the second term in ( [2.30| ), using estimate ( |2.21| ), is 

O ^ (3 33) 

uniformly in k for \k\ < /«o- 

The same method can be used to prove semiclassical asymptotics for the suscepti- 
bility x and also for higher order derivatives in k of Q. | 



Now let us denote by the energy surface at energy /i. 

J2l = {(q,p)£lR 2n --H K (q,p)=^ (2.34) 
and by da* the Liouville measure on : 

*^(9.P) = ]^| ( 2 - 35 ) 

defined for any non-critical /1 (i.e. VH K (q,p) ^ on XJD Lebesgue Euclidean 
measure). Then we have 
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Corollary 2.3 (Landau diamagnetism) Let \ be defined by ( \2.dj) , and p, be non- 
critical for H . Then for n = 2 and k = we have 



lim X = o / , II B(q) f da" (2.36) 

which is nothing but Landau's result of the diamagnetism for a 2- dimensional free 
electron gas. 

3 A "trace formula" for the magnetization 

For a temperature regime (3 > h~z +£ (e > 0), in order to get information on the 
semi-classical limit of magnetization M, we shall use the Fourier inversion formula 
instead of Cauchy formula. First of all let us remark that fp is not in the Schwartz 
space <S(R) so it is more convenient to take its derivative, which is in iS(R). We have 
explicitly 



,py ' 4cosh 2 (/3a;/2) 
So we have the following Fourier transform formula 

wt/P 



f^) = ~f dt . l( im / tX (3-1) 
1 2tt Jsr smh(7rf/p 



smh(vrt//3) 

So we can write 



where 



f'JH -p) = h- 1 I + °° dt e u ^-^ h ^ t/a / (3.2) 
' 3 sinhvrt/a v ; 



a = (3h (3.3) 

The parameter a, which has the dimension of a time, will be important in what 
follows. It plays a role in the Kubo-Martin-Schwinger condition (see [[<]]). 
Because we cannot compute the semiclassical evolution for infinite time, we shall 
consider a "smeared out" magnetization defined as follows. Fix tq and r : < tq < 
t ; given a C°° even function p such that : 



p{t) = 1 if |*| < 1 
p(t) = if |*| > 2 and / p(t)dt = 1 (3.4) 
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Let us define for r > 0, 

p T {t) = pit/r) (3.5) 



and 



M T = Tr I (f a *p T ) (^_J^\d K H K \ (3.6) 



where g denotes the inverse Fourier transform of g. Clearly M T — » M when r — > oo. 
However we shall not be able to let r — > oo in this paper (see discussion at the end 
of this section), and we will only obtain results for finite r. 

Consider a cut-off function 6 G C^(IR) with supp 6> C [—5,5], and 6> = 1 on 
[— |, |]. It is a priori arbitrary but in the sequel, we will take 5 so small that if p is 
non-critical for H K , any A G ]p — 35, p + 3<5[ will remain so. 

Let us write the following decomposition 

M T = M Tfi + M T ^ e (3.7) 

where 



' M Tfi = Tr {(/, * p r ) - /i) g> k # k } 

M ril _ e = Tr {(/, * p T ) (%^) (1 - 0) (# K - /x) <9 K # K } 



(3.8) 



We now prove 



Lemma 3.1 Let us assume (Hl-3) for H K , and let a > cr > 0. Then M Ti i_g has a 
complete asymptotic expansion in %. 

Proof . (1 — 9)(x) is supported by the union of (— oo, —5] and [5, +oo), which yields 
two contributions to M t ^_q that we call M^ = 1 _ e . 

Since f a *p T = — f£° (f' a *p T )(y)dy is the primitive vanishing at +oo of a function 
in the Schwartz class <S(R), we have for every N, uniformly for a > o" an d Ti g]0, 1], 

|Tr {(/, * p T ) (^j^J (1 " 0) + {H K - p) d K H K ^ | < C N h N 
We now consider M Tl _@. 

Since j^{f' a * p T ){y)dy = (f^ ■ p r )(0) = 1, we clearly have, uniformly for any 
A G sp{H K — p) fl (— oo, —5], and any a > a : 
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(3.9) 



and, since H K is semi-bounded from below, the contribution of 1 in ( |3.9| ) gives 
Tr |(1 — 6)~(H K — p)dH K } and obviously has a complete h expansion by the func- 
tional calculus (in fact it is of the form Tr j#i(if K — p)d K H K } for some 6\ G C^(IR).) 



□ 



The next step is to decompose p T in order to isolate the neighborhood of t = 
of the rest : 



Pr = Pro Pt + (1 - Pro) Pr = Pro + Pl,r 

since p T0 p T = p ro if r > 2r . 

This yields p T = p TQ + pi jT and correspondingly : 



(3.10) 



M Tfi = M + M o 



(3.11) 



M = Tr Uf a * p T0 ) (^-^) 9 (H K - p) d K H K 



M osc = Tr ( (f a * p liT ) ( " ^) ^ 



Finally M T is decomposed into 



M T = M + M a 



where M osc is given by ( |3.11| ), and 



(3.12) 



(3.13) 



M = M + M,i_ e 



We prove : 



(3.14) 



Lemma 3.2 Assume also p is non-critical for H K . Then for r small enough the 
classical flow induced by Hamiltonian H K on has no periodic point of period < y 
and Mq admits an asymptotic expansion in h, uniform for a = f3h > o~o > 0. 
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Proof. We have 

= - J" +3S d\ Tr Uf 9 * p T0 ) (^P^) (H K - fi) d K H K ^ + 0{h°°) 

(3.15) 

using the fact that f' a *pr is in the space <S(R) and the support property of the cut-off 
function 9. Now, if /i is non-critical for H K , it will remain true for any A e [/i, /i + 35] 
for small enough 5. 

We can rewrite the Trace inside the integral in fl3.15| ), using inverse Fourier 
transform, 

(omitting the index k of H K for simplicity). Now, using either WKB method (see 
for instance P3|), or the coherent state decomposition 0], a complete asymptotic 



expansion in H of M can be obtained, of the form : 



M = /T n+i (C (X) + ftCi(A) + ■ • • + /i fe C fc (A) + •••] mod 0(/i°°) 

which can be further integrated with respect to A on the interval [fi, fi + 35], yielding 
the result. □ 



Remark 3.3 Above Lemmas therefore imply that M has a complete asymptotic 
expansion in h. It is, so to say, the analog of the (complete h-expansion of) the 

"mean density of states" in the Gutzwiller trace formula. The other term M osc will 
be the sum of highly oscillating terms, also in complete analogy with the oscillatory 
part of Gutzwiller trace formula. Before showing this now, let us remark that the 
dominant H contribution to M has not yet been shown to reduce to the well known 

"Landau diamagnetism" . This will be postponed to the end of this section. 



Proposition 3.4 Assume (Hl-3) together with 



14 



(H.Jf.) /i is non- critical for H K (\k\ < k,q) 

(H.5) On the set (T^) T of classical periodic orbits denoted 7 with period smaller 
than r is such that the corresponding Poincare maps P 7 do not have eigenvalue 
1. 

Then for any o~\ > and for n > small enough, we have the following uniform 
asymptotics for f3Ti = o g]0, <7i] and \k\ < n , 

m osc = £ e ^/^f) ( .^ (r z\, 1/2 . r^ffi . + e + o(%°°) 

7£ ^ )t [|det(l-P 7 )|V2 smh(7rT 7 /(7) g 7 j 

(3.17) 

w/iere m 7 = J 7 dt<9 K H K (q t ,p t ). 

T* is the primitive period of orbit 7, 

5 7 (resp. i^ 7y ) is the classical action (resp. Maslov index) of orbit 7. 

f0) are constants depending on orbit 7, on the function p ljT , and on 7. 
Moreover the different orbits 7 can be choosen such that they depend smoothly on 
the parameter k and the asymptotic expansion holds uniformly in k for \k\ small 
enough. 

Proof . Since p\ )T is supported away from zero, we can rewrite M osc (defined in 
( |3.12| )) as the non-singular integral : 

Now the method that we have developed in || applies, and yields the desired result. 
□ 

Remark 3.5 At zero magnetic field (k — 0) the term m 7 is computed as follows. 
We have 

m 1 = i d K H K (q t ,pt)dt = - f adq = -$ 7 
where $ 7 zs the flux of the magnetic field through the closed curve 7. 



We now come back to M, and prove that the dominant contribution of the 
/i-expansion is indeed the well-known diamagnetic Landau term. 
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Proposition 3.6 Assume (Hl-4). Then for any o\ > and for k,q > 0, uniformly 
for (3% = a £]0, er-J we have, mod 0(h°°), 



2-n 



M=-K^I do°\\B(q)\\* + £ c k (p,0,a,T)h k + O(?i°°). (3.19) 

k>3—n 

Proof . Whereas the existence of a complete asymptotic expansion for M results 
immediately from lemmas (|3.1|) and (|3.2[ ), the explicit calculus of the dominant 
contribution in fl3.19|) is not an immediate consequence. It will be done through the 
functional calculus. In the previous section we have shown that the coefficients of 
the asymptotic ^-expansion are regular functions of k, and can be differentiated with 
respect to k. Thus, according to ( 3.14Q , we shall obtain the dominant contribution 
to M from two different contributions fio(^)U=/i an d ^ T ,i-e(X)\x=n by differentiating 
with respect to K : 



fio(A) = Tr 



(3.20) 



%,i-e{X) = Tr 



Let us consider the contribution of second derivatives in A of fto(X) : 

fi'o'(A) := G(X) = h- n c'^X) + h~ n+2 4(X) + 0(h~ n+3 ) (3.21) 

and we shall identify cq(X) and 02(A) by the following trick (inspired from ref. p3[ 
prop. V.8) : Take tp G Cq°(]// — 35, \i + 35[) and integrate against ( |3.21| ) ; we get : 



dX <p(\) G(X) = h- 1 I p eff (±) Tr 



<po[H-\ 



dX 



(3.22) 



where 



Peff(t) = Pr Q {t) 



nt /a 



sinh (nt /a) ' 



p e ff(X) is its Fourier transform, and 



ME) -=V{E) 9{E + X-fi) 



( 3.22|) follows from : 
G(X) 



dt p eff (t) Tr 



e -u(H K -x)/n e (T K -fi 



(3.23) 
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Now Q3.22I) is rewritten as 



dX <p(\) G(X) = — dX p eff (X) Tr ^ e (H K - Xh 

Z7T 



(3.24) 



which can be developed through Taylor's formula (since integration variable A is in 
a compact interval), as : 



dX ip{X) G(X) 



^ E { -^T- I d\ \ k Pe//(A)Tr \4\H K 

k=0 



oo ikh k 
k=0 



(3.25) 



Actually the term with k = is absent since p' e //(0) = (p is an even function, and 
so is p e ff)- We calculate the coefficients of K~ n and h 2 ~ n by the functional calculus, 
like in section 2 : 



dX<f(X)G(X) 



ip B (H(q,p))dp d<? - — / <p'o(H(q,p)) 



- y Pe//(o) / ^ (Ml ^ rfp + 0(/> 3 



II 5 



dq dp 
(3.26) 



Clearly the first and third terms in (|3.26| ) are independent on k (through the change 
of variable p i— »■ p — Ka(q)) , and we are left with lower order term 



h 



2-n 



12 • 4vr 2 
_^2-n 



dq dp ip'g(H(q,p)) 
d 2 

dX (pe(X)- 



k 2 II B || 2 



K II 5 



E ^ 

l<j,k<n 



E fi?, fc v)d^(g,p)(3.27) 

l<j',fc<n 



12-4vr 2 J dX 2 Jj2" x 
where we have used integration by parts : 

/j?2n G(q,p) <p"((H K (q,p))dq dp = J dX j , 

Therefore since the above calculation holds for an arbitrary test function ip, we can 
identify the functions c and c 2 (A) appearing in ( |3.21| ), modulo ^-independent terms 
as : 



da K x G(q,p) 
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f c (A) = 

2 (3.28) 
c 2 (A) = -W I ^/ ||5f^A(g,p) 

Z^A 

We can do the same calculus for fl Ty i-g instead of f2 by replacing To by r and 9 by 
1 — 9. This yields a contribution to the magnetization which, added to that coming 
from 02(A) in (|3.28|) gives the dominant Landau term in (|3.19|) . 

We shall extend now the above results to the magnetic susceptibility \- The 
statement is the following 

Theorem 3.7 Let us assume H-l to H-5 and a = (3% g]0, o~x] where o\ > is fixed. 
For Xt = X * Pt, i~ > 0, we have the decomposition 



with 

X 



h 



2-n 



24vr 2 JJ2 



Xr = X + Xosc (3.29) 



J da;\\B(q)\\ 2 + c x ^K,a,T)h k + 0{h°°) (3.30) 



k>3—n 



AS-jn+u^) \ Pi,r(T 7 ) r^/2a 



X ° SC J?)r \|det(l-P 7 )|V2 sinh(7rT » + J 

0(^°°) (3.31) 



where we have used the notations in proposition (\3.4i) and r 7 = c Xi7; d x>7 are 



7 



smooth coefficients depending on the periodic orbit 7, on a, and on function p. 

Proof We use the same cut-off already introduced for the magnetization M. So we 
define in a natural way 

X = d K M Tofi + d K M T ^ e (3.32) 

Xosc = Xt-X (3-33) 

Compute first the term Xr ,e = d K M TOi g. From the proof of proposition ( |3.4| ) we get 

1 dX r dt -^-p (±) d K Tr (e-f 9(m -X) d K % 

(3.34) 



a Jo J-00 smh(7rc/(7j \r 
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We compute derivative in the parameter k with the following easy consequence of 
the Duhamel formula 

<9 K Tr { e -^-"> 6(H-n) d K h} = Tr d K [6 (h - X) d K H]} + 

^{f Q ds (e^dj^e-^ e-^ r ^9 (H - p) d K X) (3.35) 

Then due to the support property of r , the only stationary points corresponds to 
the period T = and the leading term in fi is given by the first term. The term 
Xt,i-0 = d K M Ty i-g is computed in the same way and the both term combines to yield 
the asymptotic expansion of %. 

For the term Xosc we start from a formula like (|3.34|) replacing the time cut-off p TQ by 
the following pi )T = p T (l — p T0 ). Hence applying the methods of f| we can compute 
with the stationary phase theorem the contributions of the periodic trajectories with 
period T 7 e (IV) T . 

Remark 3.8 In the so-called "mesoscopic regime" examined in this section (i.e. 
T = for some fixed a having the dimension of time) , and in the special case 
of dimension 2, the dominant semi-classical contribution to M and Mi to M osc 
are of the same order (apart from highly oscillating factors). A comparison of the 
corresponding contributions xl and Xi t° the susceptibility is made in the physics 
literature, measuring a factor of 100 for Xi/Xl W<\h 



Remark 3.9 Thus the magnetic response is a measurable quantity where the skele- 
ton of the periodic orbits of the classical motion manifests itself clearly ; we have in- 
vestigated this effect rigorously and in great generality. Furthermore the oscillations 



in ( \3.11j ) are a generalization of the well-known de Haas-van Alphen oscillations of 
the magnetic response which are a result of the classical cyclotronic orbits demon- 
strated in dimensions 2 and 3, and which can be recovered from ( \3. 1 % ) in the limiting 
case where all classsical orbits are of cyclotronic nature (V = or quadratic). 

Now, we want to comment about the fact that we have only been able to give 
semi-classical expansions for "smeared out magnetizations" M T instead of the true 
one (r = oo). For non-zero temperature T ^ 0, the exponential decrease of /^(t) 



when k — > +oo lets us expect that the Fourier inversion formula (3T) combined with 



'trace formulas" will be enough to obtain Proposition (3~3) without the p T which 
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cuts off time at \t\ < r. We expect that our method using semi-classical evolution 
estimates for coherent states f| will allow to prove this for a = (3fi > cr > with 
suitable assumptions on the classical flow. This is presently under study. However 
for T = 0, the cut-off p T will be necessary to make the sum over periodic orbits finite 
and thus convergent, and we cannot expect to get rid of it. 

For the moment, using estimates proved in ||, we can see that it is sufficient to 
control the periods of the classical flow in the time interval [r, Co log(^]). In || we 
have proved that the semi-classical propagation of coherent states is valid in time 
interval [— c log(^), Co log(^)] for some Co > 0. So we can write down the operator 
e -f (-H-A) ag a p our j er integral operator with a complex phase for \t\ < colog(^). So 
we have to compute two terms, (H K = H), 

F^a) := J^dtTr L-^^-^ 9 (H - //) Ap 



F 2 {h, a) := J r dm je"^^ 6 (H - fj) A 



1-p 




where A is some quantum observable and R a (t) = si ^J t / a ■ The term Fi is difficult 
to check and we have nothing to say about it here except that for each time, it 
is a Fourier integral with a known complex phase but it is difficult to control the 
stationary phase argument for large times. 

The term F 2 is easily controlled because it contains the damping factor R a . More 
precisely we have 

Lemma 3.10 There exists C > such that for every % G]0, 1] and o > we have 
easily : 

F 2 (h,a)<Cc \og(±-)h^ C0 . (3.38) 
So that F 2 (h, a) is negligible for ^ large enough. 

2 

4 The regime of temperature 

In section 2 we have shown that the functional calculus applies to the ther- 
modynamical functions in the grand-canonical ensemble and provided asymptotic 

2 

expansions in the semi-classical limit provided T > h^~ £ (some e > 0). In section 
3 we have investigated a rather different temperature regime (called "mesoscopic" ) 
where ksT = h/a (some a > but finite) where a splitting of the magnetic response 
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into a "mean part" and an "oscillating part" appears in the semi-classical limit. In 
order to be complete, the "in-between regime" is now considered. 



Theorem 4.1 Assume (HI. 4)- Then the magnetisation M = d K Q has for any 

temperature T satisfying Ti l ~ £ < T < Ti^ e (some e > 0) a complete asymptotic 
expansion in h obtained by taking the derivative in k of the formal expansion in ft 



for Vt given ( 2.14 )- 



Proof . As in section 3 take r > so small that, the classical flow induced by H K 
has no periodic point with non-zero period G [— 2r ,2r ], and take p TQ as in section 
3. Futhermore let 9 G Cq°(IR) be, as in section 3 (6 = 1 on 
IR \ [—5, 5}). We decompose M and 



6 5 
21 2 



and e on 



M = M e + Ml. 



with 



(4.1) 



and similarly for M\-q. 
Furthermore : 



Tr [f p (H R - /i) 6 (H K - n) d K H K ) 



(4.2) 



Ma 



- J°° dX Tr {ft (H K - A) 6 (H K - fi) d K H K ) (4.3) 

- r dx\ r dt^I^Tr le«B-» (h _ ,) q^X 

nJ-oo sinn ixtjo I v ' J 

M e , p + M e ,i- P (4.4) 



where we insert, inside the integral over t, the partition of unity 

1 = p T0 (t) + (1 — p T0 )(t), which yields, correspondingly a splitting of M g into the two 

contributions. ■ 



Lemma 4.2 Assuming (H.l-3), Mi_g has a complete asymptotic expansion in h. 



Proof . We can proceed as in the proof of Lemma ( |3.1| ), by splitting (1 — 6){x) into 
the sum of two disjoint functions (1 — 9)^ supported respectively in [5, +oo) (for + 
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sign) and (—00, —5]. Since fp is the primitive vanishing at +00 of a function in the 
Schwartz class of C°° functions of rapid decrease, we have 



|Tr [ff> (H K - n) (1 - 9) + (H K - n) djT K ) \<C N h N (for any N) 



and 



bh' 

dn 



H - fi)^}\<C N h N 



for sany N) 



Finally, we know, like in the proof of Lemma (3J), that Tr |(1 — 9) (H — fi)d K H^ 
has a complete asymptotic fi expansion by the functional calculus. | 

Lemma 4.3 Assuming (Hl-3), one has, for Ti l ~ £ <T< %^' e 

M e ^ p = O (e-^) 
where c% is a positive constant only depending on tq. 

Proof . Using (|4.3|) , the support property of 9, and the exponential decrease of fL 
it is easy to show that : 



M e = - J^ +2S d\ Tr {/£ (H K -\)e(H K - //) dH R ) + O (e~^ 



Therefore 



rfi+28 r+00 

Me,i- P = - dXh" 1 dt(l-p T0 (t)) 

Jfi J -00 smh irt/a 

Tr L-m-xyn e ^ _ ^ Q Jj\ + ( e -8/VK\ 



nt/a 



and since, in the considered temperature regime 



id, /a 



< c e - c o\t\/h e 



sinh(7rt/(T 

we have, using the support property of 1 — p TQ : 

\M e ,i- P \ < C e-^ h " 
c\ being a positive constant depending on t . 
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Lemma 4.4 Assuming (H.l-4), then Mg tP has a complete asymptotic expansion in 

h. 

Proof . As is the previous section, we take 5 so small that, if ji is non-critical for 
H K) then any A G [fi, fi + 26} is also non-critical for H K . 

Now, using the support property of p TQ , and either WKB method, or decompo- 
sition over coherent states, a complete asymptotic expansion can be obtained for 

£> » (H, - „) a} 

for any A G [//, /x + 25]. Integrating with respect to A in this interval yields the result. 
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